Excitation of confined modes on particle arrays 
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We describe both theoretically and experimentally the existence and excitation of confined modes 
in planar arrays of gold nanodisks. Ordered 2D lattices of monodispersive nanoparticles are manufac- 
tured, embedded in a silica matrix, and exposed to evanescent prism-coupling illumination, leading 
to dark features in the reflectivity, which signal the presence of confined modes guided along the 
arrays. We find remarkable agreement between theory and experiment in the frequency-momentum 
dispersion of the resonances. Direct excitation of these modes reveals long propagation distances 
and deep extinction features. This characterization of guided modes shows the great potential of 
metallic particle arrays for optical signal processing and distant sensing applications. 

Metal nanoparticles are capable of producing high levels of confinement and enhancement of incident light pQ. 
These appealing properties make them useful for sensing (e.g., through detection of plasmon-resonance shifts [2HS] 
and surface-enhanced Raman scattering [7HTT]). cancer therapy [12J[T3], drug delivery [14], improved photovoltaics 
[15] [16] , and catalysis [17] [18] . These applications have attracted a considerable degree of interest on the optical 
properties of nanoparticles, which strongly depend on size, morphology, and composition [19-21 ]. The interaction 
between neighboring particles can also produce changes in their optical behavior, which are dramatic when small 
particle gaps are involved [22 -24]. Additionally, long-range interactions lead to collective modes in ordered [25) [26] 
and disordered [27] particle arrangements, which have been proposed for waveguiding in both ID [28] [29] and 2D 
[25] [30] configurations. 

In particular, lattice resonances in 2D ordered particle arrays act as an effective way of creating guided modes 
[3T] that are strongly dependent on the choice of materials, particle size and shape, and lattice parameters, leading 
to distinct TE and TM modes [25] [32] . This configures a versatile platform for applications to distant sensing and 
waveguiding, which deserve a detailed characterization of the arrays. Direct far-field sampling of quasi-bounded modes 
in 2D arrays through diffracted orders have been extensively reported [33-38 . Additionally, theoretical descriptions are 
available [39], which become particularly intuitive when assuming dipole inter-particle interactions [25] [31]. However, 
guided modes in metal particle arrays have not been yet directly sampled (i.e., sampled without the assistance of 
diffraction to couple to far- field radiation), and a detailed experimental study of their optical behavior is still missing. 




FIG. 1: Images and scheme of the system under consideration, (a) and (b) SEM images of our nanoimprinted arrays. The gold 
nanodisks are arranged in a square lattice with parameter a — 300 nm. The height of the disks is h — 50 nm and the diameter 
is D — 200 nm in (a) and D — 250 nm in (b). (c) An array of golden disks is embedded in a homogeneous SiCb environment 
(refraction index nsio 2 — 1-46) at a distance d — 600 nm from a prism interface. This is the Otto configuration, whereby light 
undergoes total internal reflection inside the prism (uf 2 = 1-67) with incident angle 0i and parallel wave vector k\\ outside 
the light cone. Light couples evanescently to the confined modes of the array, thus producing a measurable attenuated total 
reflection. 
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FIG. 2: Simulations of confined modes in a square array of SiCb-embedded gold nanodisks of D = 200 nm in diameter, (a) 
Dispersion diagram showing the light energy and wave- vector dependence of the array reflection coefficient (color scale). Guided 
modes show up as bright features. The light cone (solid line) and its first Bragg- reflect ion (dashed line) are plotted for reference, 
(b) Electric- field intensity maps over a unit cell at a distance D/2 above the array for characteristic points within the dispersion 
diagram, as shown by labels A-D in both (a) and (b). The nanoparticle contour is represented by a dashed circle. Dipolar 
patterns are visible when confined modes are excited (A and B), leading to large field enhancement (see color scale). In contrast, 
lower field enhancement is associated with non-resonant reflection (C and D), although hotspots can localize between particle 
gaps (see D), leading to an almost dispersionless feature in (a). 

In this Letter, we report direct observation of guided modes in ordered gold-particle arrays through dips in their 
evanescent-field reflection using the Otto configuration [40j ST]. Gold nanodisks have been fabricated by imprint 
lithography [42] and embedded in a symmetric silica environment [Fig. [TJa,b)]. A high- index prism is then used for 
evanescent coupling of incident light [Fig. [TJc)]. The reflected signal unambiguously shows the presence of guided 
modes. We further report electromagnetic simulations for both the guided-mode dispersion relations and the reflec- 
tivity under prism illumination. Experiment and theory are in excellent agreement for the dispersion relations and 
the mode propagation distances. Our results bear some resemblance to guided modes in photonic crystal waveguides 
[43U45] and electromagnetic surface waves in drilled metal films [46] . 

The array modes under study originate in the lattice resonances produced by long-range interaction between the 
particles close to the onset of different diffraction orders [25 j. Under such conditions, there is constructive interference 
of the fields scattered by the particles, leading to array-confined modes, which are in fact guided modes that propagate 
along the plane of the array. These modes show up at frequencies to the red with respect to the light line in a frequency- 
wavevector diagram 

[31]. In Fig.ga) we show the calculated reflection coefficient of arrays of gold nanodisks placed 
in a symmetric Si02 environment. The dispersion diagram is dominated by two types of features: a prominent one 
arising from the noted lattice resonances; and a weaker one associated with the excitation of localized surface plasmons 
that are confined at the gaps between neighboring nanoparticles [23] (see below). The reflection coefficient is peaked 
along a line close to the light cone (solid line), indicating the position of the dispersion relation of the lattice modes. 
For reference, we also show in the diagram the umklapp folding of the light cone (dashed line), which crosses the light 
cone at the boundary of the first Brillouin zone (BZ). While the lattice resonance within the first BZ is both too 
narrow and close to the light line to be resolved in the figure, it is clearly visible withing the second BZ. 

The intensity of the electric field, which is represented in Fig. [2^b) for a unit cell at a distance D/2 from the 
middle plane of the lattice, clearly shows the collective character of the resonance. Along the dispersion relation of 
the propagating mode, such at positions A and B, we have large values of the field enhancement spreading all over 
the plane, with relatively small spatial variations. We also represent the field at the weaker, relatively fiat spectral 
feature below the resonance at point C. Clearly, the latter shows strong localization at the particle gaps, which is 
consistent with an interpretation in terms of localized gap plasmons between neighboring nanoparticles and explains 
the observed vanishing group velocity. Under non-resonant incidence conditions (i.e., away from the propagating 
mode and localized plasmon regions), the enhancement is featureless and weak, as shown for point D. 

Arrays embedded in a symmetric environment thus display sharp resonances that can be easily modeled theoretically, 
but their fabrication and characterization is difficult compared to arrays supported on a substrate. A high degree of 
asymmetry is known to strongly modify and even prevent the existence of the lattice resonances [25j [39] . In order 
to reach a compromise between these two extreme situations, we have fabricated gold nanodisk arrays on a Si02 
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FIG. 3: Experimental (a,b) and theoretical (c,d) energy-angle-dependent reflectance of the combined array-prism system [Fig. 
[TJa)] for nanodisks of diameter D — 200 nm (a,c) and D — 250 nm (b,d). The excitation of a guided mode is revealed by a dark 
(low reflectance) feature at the right of the critical angle (0i « 64 deg) . 

substrate and added a dielectric capping of the same material Qa)). This capping of thickness d = 600 nm not only 
ensures the existence of the resonant modes, but also provides protection for the nanostructure in potential device 
applications. We have measured the reflectivity spectra of these samples in a rotation stage using a coupling prism of 
F2 glass. When light is reflected at the prism-sample interface under total-internal reflection conditions, it produces 
an evanescent field that, for values of the parallel momentum matching those of the guided modes, is coupled into the 
latter. This coupling transfers energy that is dissipated due to losses in the metallic nanoparticles, therefore producing 
a dip in the reflection spectrum, which is a signature of the excitation of a surface mode [40] . 

The spectral dependence of this dip can be clearly resolved in Fig. [3J where it shows up to the right of the critical 
angle (0 C ~ 64). The agreement between theory [Fig. [3^c,d)] and experiment [Fig. [3ja,b)] is excellent in position, 
intensity, and width of the spectral features for samples with two different particle diameters [D = 200 nm in Fig. 
[3ja,c) and D = 250 nm in Fig. [3]^b,d)]. Some discrepancies are observed in the features at the left end of the angular 
region under investigation, which we attribute to systematic errors in the characterization setup, which is optimized 
for data collection to the right of the critical angle. As expected from the calculations of Fig. [2] for a symmetric 
environment, the resonance appears close to the light line and nearly parallel to it, specially for the smallest particles. 
Incidentally, almost perfect coupling takes place, with the reflectance reaching values close to zero. 

A more detailed characterization of the modes is presented in Fig. |4j Panel (a) shows excellent agreement between 
the theoretical (dashed curves) and experimental (solid curves) dispersion relations for both particle sizes. The degree 
of confinement, indicated by the separation from the light (dotted) line in Fig. [4ja) and also by the penetration of the 
field intensity outside the array in Fig. |4^b), is also in excellent agreement. Notice that this penetration distance [37] 
is comparable to the separation of 600 nm between the prism and the particles, but the modes are still well defined 
and clearly resolvable. 

Finally, we have calculated the propagation distance along the plane of the array [Fig. |4jc)] by fitting the re- 
flectance dips to Lorentzians [see insets in Fig. [4^a)]. Theory predicts larger propagation than what is observed in 
the experiment, and we attribute this discrepancy to imperfections in the shape and position of the particles. Actu- 
ally, the discrepancy between theory and experiment is less pronounced for the largest particles under consideration 
(gray curves), for which variations in morphology with the same degree of absolute tolerance are relatively smaller. 
Propagation distances of a few microns (~ 10 surface-mode wavelengths) are obtained. 

In conclusion, we have experimentally demonstrated the existence of confined propagating modes in planar arrays 
of plasmonic nanoparticles, and we have determined their dispersion relation and characteristic lengths. Our electro- 
magnetic simulations are in excellent agreement with the experiments. These results should trigger further research 
into the direct spatial mapping of guided modes. Our particle arrays can be potentially used to propagate optical 
signals. They can also assist distant sensing, whereby the optical signal from the object or substance to be sensed 
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FIG. 4: Experiment-theory comparison of the dispersion relation (a), the penetration length outside the array (L z ) (b), and 
the propagation distance (L x ) (c) for the guided modes obtained from the position and width of a Lorentzian function fitted 



to the data of Fig. [3 
reflection dip in Fig 



Theory: dashed curves. Experiment: solid curves, obtained from constant-energy Lorentzian fits of the 
3] [the inset in (a) shows a characteristic example of the Lorentzian fit for the energy 1.74 eV signaled with 



crosses in the main plot] . Black (gray) curves correspond to disks of diameter D = 200 nm (D = 250 nm) 



propagates through an array to the position at which the detector is placed. This possibility provides a natural hybrid 
between plasmon-based sensors and nanoparticle-based waveguiding. 
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